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We describe our explicit Lorentz-invariant solution of the Einstein and null geodesic equations for 
the deflection experiment of 2002 September 8 when a massive moving body, Jupiter, passed within 
3.7' of a line-of-sight to a distant quasar. We develop a general relativistic framework which shows 
that our measurement of the retarded position of a moving light-ray deflecting body (Jupiter) by 
making use of the gravitational time delay of quasar's radio wave is equivalent to comparison of the 
relativistic laws of the Lorentz transformation for gravity and light. Because, according to Einstein, 
the Lorentz transformation of gravity field variables must depend on a fundamental speed c, its 
measurement through the retarded position of Jupiter in the gravitational time delay allows us to 
study the causal nature of gravity and to set an upper limit on the speed of propagation of gravity in 
the near zone of the solar system as contrasted to the speed of the radio waves. In particular, the v/c 
term beyond of the standard Einstein's deflection, which we measured to 20% accuracy, is associated 
with the aberration of the null direction of the the gravity force ( "aberration of gravity" ) caused by 
the Lorentz transformation of the Christoffel symbols from the static frame of Jupiter to the moving 
frame of observer. General relativistic formulation of the experiment identifies the aberration of 
gravity with the retardation of gravity because the speed of gravitational waves in Einstein's theory 
is equal to the speed of propagation of the gravity force. We discuss the misconceptions which 
have inhibited the acceptance of this interpretation of the experiment. We also comment on other 
interpretations of this experiment by Asada, Will, Samuel, Pascual-Sanchez, and Carlip and show 
that their "speed of light" interpretations confuse the Lorentz transformation for gravity with that 
for light, and the fundamental speed of gravity with the physical speed of light from the quasar. For 
this reason, the "speed of light" interpretations are not entirely consistent with a retarded Lienard- 
Wiechert solution of the Einstein equations, and do not properly incorporate how the phase of the 
radio waves from the quasar is perturbed by the retarded gravitational field of Jupiter. Although all 
of the formulations predict the same deflection to the order of v/c, our formulation shows that the 
underlying cause of this deflection term is associated with the aberration of gravity and not of light, 
and that the interpretations predict different deflections at higher orders of v/c beyond the Shapiro 
delay, thus, making their measurement highly desirable for deeper testing of general relativity in 
future astrometric experiments like Gala, SIM, and SKA. 
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1. INTRODUCTION 



The recent Very Long Baseline Interferometric (VLBI) experiment of the uhra-precise measurement of the relativis- 
tic deflection of a quasar's radio waves (light) by the gravitational field of Jupiter has captivated the keen interest 
of relativists who were trying to deeply understand gravitational physics. The experiment and its physical interpre- 
tation has illustrated some of the important subtleties in the understanding, the approaches, and the mathematical 
techniques of general relativity (GR). 

The experimental goal was the measurement of the (non-radial or tangential) deflection component associated with 
the retarded position of Jupiter caused by its orbital motion and mathematically connected to the finite value of 
the fundamental speed entering in front of the time derivatives from the metric tensor in the gravity field and null 
geodesic equations (that is, the fundamental constant c in the Christoffcl symbols and the curvature tensor). This 
fundamental speed of gravity causes the retardation in propagation of gravity force from a moving massive object 
which can be measured in the light-ray deflection experiments. In this paper we use general theory of relativity and 
denote the fundamental speed of gravity by c. Label Cg for the speed of gravity is understood as a parameter in 
the experimental fitting of observations to the Einstein theory. The symbol c is frequently called as the speed of 
light. However, one must clearly separate the relativistic effects caused by the physical speed of light from those 
associated with the fundamental speed of gravity in the Eiunstein equations. We shall show in the present paper how 
to distinguish these effects to avoid ambiguity. 

Einstein's general principle of relativity assumes that the gravity field obeys to the special principle of relativity 
and transforms from one inertial frame to another in accordance with the Lorentz transformation of the Minkowski 
space-time. Hence, it postulates that in a geometrized system of units '3| the speed of gravity c = 1 in any frame 
of reference. We measured the impact of the retarded gravitational field of Jupiter on light deflection (as well as the 
instantaneous Shapiro deflection of 1100 /larcsec) of 50 ± 9 /larcsec, in excellent agreement with the GR prediction 
of 51 /iarcsec based on assumption that c — 1 irrespectively of the reference frame. We have interpreted these result 
as a measurement of the fundamental speed of gravity c = 1.06 ± 0.21 and confirmation that gravity interacts with 
moving particles (photons) not instantaneously. 

However, at present there are four other interpretations, which do not agree among themselves, concerning the 
results of the experiment ^ . 

• the experiment measured the speed of the radio waves from the quasar in the static gravitational field of Jupiter 

• the experiment measured the aberration of the radio waves from the quasar |^, Q , 

• the experiment measured no significant property of light or of gravity |^ , 

• the experiment measured the light-time delay discovered by R0mer in 1676 

In this paper, we streamline our previous theoretical calculations by working out an explicit Lorentz-invariant 
Lienard-Wiechert solution of the Einstein gravity field equations in §2. The VLBI measurement of the phase of a 
radio wave, which has propagated in the field of a moving body (Jupiter), is described in §3. We demonstrate that 
the fundamental speed of the characteristics of the gravity field equations (which we call for brevity "the speed of 
gravity") and the aberration of gravity force are indeed associated with the experiment. The major misconceptions 
about the experiment and properties of gravity in general relativity are summarized in §4. Finally, in §5 we outline 
why these alternative interpretations of the experimental results differ with our general relativistic interpretation, even 
though all interpretations predict the same angular deflection result to order v /c beyond the Shapiro delay. Appendix 
describes how light propagates in a simple bi-metric model of gravity with two metrics, and provides a clear evidence 
that a light-ray deflecting body deflects light from its retarded position connected to observer by a null line lying on 
the gravity null cone that is because gravity propagates with finite speed. 



^ Recent publication of Frittelli |^ agrees witii our understanding of the gravitational physics of the Jovian deflection experiment 
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2. THE METRIC TENSOR, RETARDED TIME AND THE FUNDAMENTAL SPEED 

2.1. Formulation 

Our interpretation of the experiment is based solely on General Relativity and the physical meaning of the retarded 
Lienard-Wiechert potentials used to solve the Einstein equations. These potentials describe the propagation of the 
gravity field even in the case of gravitating bodies which move uniformly with constant speed [I^ exactly as in 
electro dyna mics where the retarded Lienard-Wiechert potentials describe electromagnetic field of a moving point 
charge fill IT^ . We emphasize that although gravitational waves, which amplitude is falling off as 1/r, are not 
generated by a uniformly moving body the null characteristics of the gravitational field are precisely defined by the 
Lienard-Wiechert potentials through their property of the Lorentz invariance and the principle of physical causality. 
Thus, our use of the "speed of gravity" term is a synonym for the ultimate speed of the null characteristics of the 
gravity field equations. This speed is numerically equal to the fundamental speed of the Minkowski geometry which 
is associated with the speed of light in vacuum. However, in general relativity the fundamental speed could be 
diff'erent from the speed of light. In the near-field zone of a gravitating system the null characteristics of gravity reveal 
themselves through the Lorentz invariance of the gravitational force (the aberration of gravity) when one compares 
observations done in one inertial frame with another, while in the far-field radiative zone the null characteristics of 
gravity can be traced as freely propagating gravitational waves. 

In general relativity light and gravity null rays are bi-characteristics, that is they propagate with the same funda- 
mental speed. At the first glance it may look that they are identical and can not be discriminated. However, this 
point of view is erroneous. Although light and gravity propagate with the same speed on the hyper-surface of one 
and the same null cone, they (1) generally propagate in different directions in space in each particular gravitational 
experiment, and (2) they are associated with different physical effects (the aberration of light and the aberration of 
gravity, for instance). Hence, discrimination between the light and gravity null cones is possible if one makes use of 
additional properties of the propagating system of the light and gravity null rays. 

Our definition of the "speed of gravity" is more general than that used by Asada, Samuel, and Will IE IS 
who limited its meaning as the speed of propagation of free gravitational waves alone. In their formulations of the 
experiment, these authors assumed only far-field gravitational effects, where gravitational waves are dominant and, 
hence, the differentiation between the speed of gravity c and that of light occurs only at orders of {v/cf' beyond the 
Shapiro delay and higher. This was one reason why the "speed of light" was interpreted as causing the observed 
aberration of gravity force of Jupiter. However, the experiment was performed in the near-field of the quasar radio 
wave- Jupiter interaction where the transformational degrees of freedom for gravity not associated with gravitational 
waves are dominant but they must not be confused with the transformational degrees of freedom for radio waves used 
for observations. 

It is clear that there are different theoretical frameworks to interpret the Jovian deflection experiment as a new 
test of non-stationary properties of the gravitational field in the solar system. However, our point of view is that 
the most natural interpretation should rest solidly on GR which has passed all other experimental tests in strong 
compliance with our current theoretical understanding of gravitational physics. In other words, in order to interpret 
the experiment properly, the space-time properties which GR postulates must not be violated; that space-time is a 
differentiable manifold endowed with an affine connection (Christoffel symbols) whose geodesies form a privileged set 
of world-lines in space-time, and their knowledge (observation) allows us to extract information about the curvature 
tensor (tidal gravitational force) without ambiguity. The Christoffel symbols define the law of motion of test particles 
and extended bodies and, hence, generalize the concept of the Newtonian gravitational force in GR. They arc formed 
from partial derivatives of the metric tensor that defines geometric properties of the space-time and the causal structure 
of the null cone. At the same time the metric tensor is associated with a gravitational potential whose properties are 
determined by Einstein's equations. 

In what follows the Greek indices a,/3, ... run from to 3, the Roman indices ... run from 1 to 3, repeated 
Greek indices mean Einstein's summation from to 3, and bold letters a = {a} ,a? ,a?),b — {b^,b^,b^), etc. denote 
spatial (3-dimensional) vectors. A dot between two spatial vectors, for example a ■ b = a^b^ + a^b^ + a'^b^, means the 
Euclidean dot product, and the cross between two vectors, for example a x b, means the Euclidean cross product. 
We also use a shorthand notation for partial derivatives da — d/dx". Greek indices are raised and lowered with full 
metric gafj. The Minkowski (fiat) space-time metric rjap — diag(— 1, -1-1, -1-1, -1-1). 

The solution of the Einstein equations, using the Lienard-Wiechert potentials has been essentially understood 
and developed by previous researchers and we do not intend to go into complete details of this iterative procedure 
called by the post-Minkowskian approximations [T^ . For our purposes the first (linearized gravity) post-Minkowskian 
approximation is sufficient since the gravitational light-ray defiection experiments are not yet sensitive enough to 
measure non-linearities of the Einstein equations. The linearized gravity approximation is as follows. 
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We introduce the weak-field decomposition of the metric tensor 

ga/l = r]af3 + hap , (2.1) 

where hap is the post-Minkowskian perturbation of the Minkowski metric tensor "qap- We impose the harmonic gauge 
condition on the metric tensor 

a^/i"/3 _ = . (2.2) 

The linearized Einstein equations in the first post-Minkowski approximation are of the hyperbolic-type (wave) equa- 
tions for the metric perturbation. In arbitrary harmonic coordinates a;" — {ct,x), the Einstein equations read ^ 

Here (□ = —cr^d^/dt^ + V^) is the wave operator in flat space-time, and T^'^ is the stress-energy tensor of point-like 
light-ray deflecting bodies in linearized approximation 

N 

r^''(t,a;) =^A47-i(i)<(i)<(t)<5(3)(a;-a;Ji)) , (2.4) 

a=l 

where the index a = 1, 2, TV enumerates gravitating bodies of the solar system, Ma is the (constant) rest mass of the 
ath body, Xa{t) is time-dependent spatial coordinate of the ath body, Va{t) = dXa{t)/dt is velocity of the ath body, 

= 7a(l, Va/c) is the four-velocity of the ath body, 7^ = (l — v^jt?) is the Lorentz-factor, and i5('^)(a;) is the 
3-dimensional Dirac's delta- function. We have neglected \/ g in Eq. 1)2. 4|l because in the linearized approximation 
\J—g = 1 -l- OiK) and the terms proportional to the metric tensor perturbation are irrelevant in T^^ . 

Because the Einstein equations (|2.3|l are linear, we can consider their solution as a linear superposition of the 
solutions for each body. It allows us to focus on the relativistic effects caused by one body (Jupiter) only, although in 
the actual experiment the gravitational fields of Sun and Earth are also important and were included in the analysis 
(see ^, 14] for more detail) . Solving Einstein's equations (|2.3|l by making use of the retarded Lienard-Wiechert tensor 
potentials [l5|, one obtains the post-Minkowski metric tensor perturbation for Jupiter (isl IT^ 

/.^-(^,.)^^^f + (2.5) 

where = —Uar°', = x" — Xj{s). The metric perturbations for Sun and Earth look similar to Eq. H2.5|l with a 
corresponding replacement of masses, coordinates and velocities for the Sun and Earth respectively. In Eq. (|2.5(l . Mj 
is the mass of Jupiter, = c~^dx'^{s)/dT = j{s){l,vj{s)/c) is its four-velocity, with r = r(s) being a proper time 
on Jupiter's world line, vj{s) = dx{s)/ds is Jupiter's coordinate velocity, and 7(5) — ds/dr = (1 — Wj(s)/c^)~^/^. 
Notice that the metric tensor perturbation (|2.5|l is valid for accelerated motion of Jupiter and is not restricted by the 
approximation of a body moving on a straight line (see [Tsf for more detail). In other words, the four- velocity u" in 
Eq. 1)2. 5|l is not a constant. 

Because we solved the Einstein equations 1)2. 3|l in terms of the retarded Lienard-Wiechert potentials, the distance 
r°' — — Xj{s)jthe Jupiter's worldhne Xj{s) = {cs,xj{s)), and the four-velocity u"(s) are all functions of the 
retarded time s [15j. The retarded time s is found in the first post-Minkowski approximation as a solution of the 
gravity null cone equation 

Ty^.r'^r^ ^ r^^, (x^ - :r^(s)) (.t^ - ^[^(s)) = , (2.6) 

that is 

s = t--\x-xj{s)\, (2.7) 
c 

where the fundamental constant c in Eq. (|2.7|) must be conceptually understood as the speed of characteristics of 
the gravity field equations because it comes from the wave operator of the Einstein gravity field Eqs. (|2.3(l . Indeed, 



^ t is time and x is space coordinates. 
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Maxwell equations have not yet been considered; hence the retardation in Eq. (|2.7() cannot be associated with the 
propagation of electromagnetic waves (light). In what follows, the retarded time s will appear in all subsequent 
equations related to the propagation of radio waves. It is important to keep in mind that Eq. (|2.7|l connects the 
point of observation x and the retarded position of Jupiter x{s) by a null line which is always a characteristic of the 
Einstein field equations. Radio waves arc also propagating along a null direction connecting the observer and the 
source of radio waves (quasar). However, the null characteristic of the Einstein equations (|2.7|l is well separated in the 
sky from the null direction associated with the propagation of the radio wave so that they can not be observationally 
confused (see Appendix A for additional details). 



2.2. Interpretation of the Einstein equations 

The retardation shown in Eqs. 12.5|l - 12. 7|) is physically associated with the finite speed of propagation of the null 
characteristics of the gravity field. In the most general case, this propagation is realized in the form of a traveling 
gravitational field which can be decomposed into components of different algebraic types |l7j . In the wave zone of the 
solar system, the traveling gravitational field decouples from the system and propagates as a transverse-traceless (TT) 
gravitational wave with amplitude decaying as 1/r where r is the distance from the system. However, in the near 
and transition zones of the solar system, a propagating gravitational field has components which are not completely 
decoupled from the matter and decay faster than 1/r (17.] . These components of the propagating gravitational field 
are closely associated with the residual gauge freedom of the metric tensor ga/3 and for this reason we call them the 
gauge modes. They are also known as transverse (T) and longitudinal (L) wave modes These gauge modes control 
the validity of the laws of conservation of matter and gravity field in the lower order post-Newtonian approximations 

MM 

All components of the propagating gravitational field of the solar system interact with radio waves moving from 
the quasar to the Earth, and perturb the phase of an electromagnetic waves according to the retardation Eq. 12. 7() 
which incorporates both the near-field and far-field effects of gravity as explained in §3. The use of the retarded 
Lienard-Wiechert potentials, rather than the advanced potentials, is consistent with the principle of causality, and 
the observation of the orbital decay of binary pulsars caused by emission of gravitational radiation . 

The speed c is a fundamental speed of the special theory of relativity and Lorentz-invariance of flat (Minkowski) 
space-time metric rjap. According to our definition the speed of gravity is a fundamental constant of the curved space- 
time manifold which is associated with temporal variation of gravitational field and appears as a physical quantity 
any time when one takes a derivative of the metric tensor with respect to time coordinate pH . In GR the speed of 
gravity is equal to c because at each point the curved space-time manifold has a tangent space-time which coincides 
with Minkowski space. But this postulate must be tested experimentally [20l[2ll| . 

The measured Jovian deflection is most sensitive to the retardation of the gravitational fleld caused by its Lorentz 
invariant properties. Thus, what we measured was not the speed of TT gravitational waves but the speed of gravity in 
the near-fleld zone. Since the speed of TT gravitational waves can not exceed the speed of gravity in the near-zone for 
any viable theory of gravity f23 | , our VLBI experiment sets an upper limit on the speed of propagation of gravitational 
waves as well. 



2.3. The gravito-magnetic field 

Einstein's equations in the frame attached to Jupiter do not contain time derivatives of the metric tensor because the 
Jovian field is static; hence the retardation of gravity (associated with time derivatives of the metric tensor and its null 
characteristics) cannot be measured in this frame. This condition has led many to infer that the null characteristics 
and retardation of the gravity field must remain unmeasurable after a Lorentz transformation to another frame, for 
example to the surface of the earth or the barycenter of the solar system. Thus, it is claimed, that the defiection- 
of- light measurement cannot be associated with the retardation of gravity. But, the Lorentz transformation of the 
Einstein equations l|2.3|l from static to a moving frame induces time derivatives of the metric tensor perturbations 
which can be measured in the light defiection experiments from the retardation of gravity fleld in order to determine 
the fundamental speed of its null characteristics. 

Our interpretation does not mean that the laws of gravitational physics change under a Lorentz transformation, 
but only that some property of gravity can become measurable after transformation from a static to a non-static 
frame. A similar example comes from electrodynamics in the case of a single electron which has a static electric field 
in its own frame. A Lorentz transformation to a moving frame induces a new field - the magnetic field, which was 
not observable in the static frame. Maxwell's equations, also, show that the strength of the magnetic fleld in the 
non-static frame depends on a time derivative of the electric field coupled with the speed of light. 
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The same analogy is associated with the Jovian deflection experiment. In Jupiter's frame there is only a static 
"electric-type" gravitational field. With a Lorentz transformation to a moving frame, a "magnetic-type" field, called 
a gravito-magnetic field, is induced. This field, according to the theory of general relativity, arises from moving 
matter (matter currents) just as an ordinary magnetic field arises from moving charges (electrical currents). The 
analogy is so apt that the equations describing this "magnetic-type" component of gravity can essentially be adapted 
from Maxwell's equations for electromagnetism by replacing the charge density with the mass density and the charge 
current with the mass current, with changes of numerical factors reflecting different helicities of electromagnetic (spin 
1) and gravitational (spin 2) fields 23]. Thus, in a moving frame, the translational motion of Jupiter produces the 
gravito-magnetic field which deflects the light tangentially by dragging it to the direction of motion of Jupiter [2ll[23 |. 
We can measure this tangential gravito-magnetic dragging of light and express its magnitude in terms of the speed of 
characteristics of the gravity field equations, in the analogous way that we can express the magnitude of the magnetic 
field interaction in terms of the speed of the characteristics of electromagnetic field equations [2ll |2^ . 

However, the predicted light deflection must be invariant to any Lorentz transformation. This is clearly the case 
because all of the calculations of the deflection experiment, regardless of the reference frame, predict the same 
deflection! But, the gravito-magnetic property of gravity that is sampled by the experiment does depend on the 
reference frame. In Jupiter's frame, all of the deflection is associated with the static gravitational field - but we do 
not observe in this frame. In the frame of the experiment (the barycentric of the solar system), the deflection is 
associated with a (radial) static gravitational field of Jupiter plus a small (tangential) gravito-magnetic field, due to 
its orbital velocity, which is coupled with the speed of gravity, and reveals itself through the retardation of gravity 
already in terms of order v/c beyond the Shapiro delay. More precise measurements of the deflection by Jupiter, 
which are directly sensitive to TT gravitational waves, will require a more exact solution of the Einstein equations 
taking into account Jupiter's orbital acceleration. 



2.4. The speed of gravity and the Newtonian limit 

The Principle of Correspondence of GR with the Newtonian theory (closely related to the Principle of Equivalence) 
does not provide any information on the speed of gravity because the speed of gravity is associated with the temporal 
changes of the gravitational field while the Newtonian limit is stationary. In principle, the speed of gravity should 
appear in the right side of the Einstein equations as a constant characterizing coupling between space-time and matter 
variables |23|. Thus, it seems natural to state that the speed of gravity can be already derived from the validity of the 
Principle of Correspondence between Einstein's theory of general relativity and the Newtonian theory of gravity. But, 
the metric tensor perturbation in the Newtonian limit of general relativity depends on a coupling constant k, = G/c^ , 
where G is the observed value of the universal gravitational constant. If the speed of gravity were different from c the 
difference would be absorbed by re-definition of the observed value of G in such a way that the coupling constant k 
is kept fixed 



3. THE DEFLECTION EXPERIMENT 
3.1. The electromagnetic phase 

Let us assume that the speed of light is equal to the speed of gravity. It simplifies calculations but the observable 
effects of the light and gravity propagation can be still clearly separated as we discuss later. For the doubting reader 
appendix of this paper gives more general covariant description of propagation of light rays in a bi-metric theory of 
gravity, thus, helping one to formally discern null characteristics of light from those for gravity given by Eq. H2.7|l . 

Very Long Baseline Interferometry (VLBI) measures the phase Lp of the wave front coming from a radio source. 
The phase is an invariant scalar function with respect to coordinate transformations, and it is determined in the 
approximation of geometric optics from the eikonal equation 0, 12H 

= , (3.8) 

where g>^^ = 77^'' — /i'"'. The eikonal Eq. H3.8|) is a direct consequence of Maxwell's equations 0,11^1^ and its solution 
describes a wave front of an electromagnetic wave propagating in curved space-time defined by the metric tensor in 
Eqs. (|2.1I) . H2.5|l . We emphasize that the electromagnetic wave in Eq. (|3.8|l is a test field that has no back-action on 
the properties of the metric tensor perturbation /i^^, and does not change the curvature of the space-time. Hence, the 
structure of Eq. H3.8|l is not related to the specific properties of the electromagnetic field, but is completely determined 
by the structure of space-time, or, more accurately, its geometry defined by Eqs. H2.1fl . H2.5|l ,26.] . 
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Let us introduce a co-vector Ka = da(p- Let A be an afRne parameter along a light ray orthogonal to the electro- 
magnetic wave front cp. The vector = dx°'/dX = g'^^dpf is tangent to the light ray. Eq. 1)3. 8|) expresses a simple 
fact that the vector K"" is null, that is g^j^K'^K'^ = 0. Thus, the light rays are null geodesies defined in the linearized 
post-Minkowskian approximation by equation |25| 

^ = ^d^K.^K'^K- , (3.9) 

which describes propagation of a light (radio waves) ray from the quasar to the observer (VLBI station). Wave Eq. 
(|3.8|) and ray Eq. H3.9|) are supplementary to each other and have equivalent physical content. 

A more straight-forward solution can be obtained to Eq. (|3.8() . If Jupiter's acceleration is neglected, a plane- wave 
solution of this equation is 

•fix") = fo + ko.x'^ + (kcu") In (-fc,r") , (3.10) 

where ipo is a constant of integration, fc" — c^^{v, vk) is a constant wave vector of the unperturbed radio wave such 
that rj^iiik^k'^ = 0, is its unperturbed constant frequency. The speed c in the expression for the wave vector fc" 
must be conceptually understood as the speed of light, and the unit vector k defines the unperturbed direction of 
propagation of the radio wave from the quasar. One can easily prove that Eq. H3.10|l is a particular solution of the 
electromagnetic eikonal Eq. H3.8|l . Indeed, observing that 

9ar'^ = - —do.s , (3.11) 

7 

where 7 = (1 — Wj/c^)^^/^, one obtains from the gravity null cone Eq. 1)2. 6|l 

a„s = -7— . (3.12) 

rR 

Differentiation of Eq. H3.10|l using Eqs. (|3.11|l and H3.12|l shows that the eikonal Eq. (|3.8() is satisfied. 

Eq. (|3.1Q|1 for the electromagnetic phase is clearly Lorentz-invariant and valid in an arbitrary coordinate system. 
The gravity field (transverse-traceless and gauge modes) propagates from Jupiter along the hypersurface of gravity 
null cone (|2.7|) and perturbs the phase front of the radio wave at a field point a;" (see figure ^ . Notice that the 
procedure of solving of the eikonal Eq. 1)3. 8|l docs not change the physical meaning of the retarded time s in the 
coordinate of Jupiter x°'j{s), affected by the propagation of the gravity. Eq. I|2.7|l does not depend on the wave vector 
k°' which alone contains the information about the speed and direction of propagation of the quasar's radio wave. In 
other words, a gravitating body (Jupiter) always perturbs the electromagnetic phase and deflects the light by acting 
from its retarded position xj{s) regardless of the direction of motion of the incoming photon and/or the magnitude 
of its impact parameter with respect to the light-ray deflecting body. Hence, the speed of the physical light (radio 
waves) used for VLBI observations can not enter in any way to the gravity null cone Eq. 1)2. 7|l . The light is used as 
a test particle in order to measure the retardation of gravity effect through the retarded position of Jupiter xj{s) 
which is determined in the experiment from the magnitude and direction of the gravitational deflection of light and is 
compared with the JPL ephemeris position of Jupiter obtained independently from its direct observation in a series 
of radio-tracking and/or VLBI measurements of spacecrafts which orbited Jupiter (see Fig. Ofor comparison of 
the two positions of Jupiter). 

This independence of the retarded time s from the electromagnetic wave vector fc" clearly shows that the retardation 
is not due to the propagation of the radio wave from the quasar. The radio propagation can not be responsible for 
the retardation effect in Jupiter's position because the radio wave propagates in a different direction (see Fig. [2J|. 
The retardation in Jupiter's position is a consequence of the finite value of the fundamental speed of propagation of 
gravity and exists independently of solving the eikonal equation; hence the solution of the eikonal equation can not 
change this property of the propagating gravity field. This conclusion is supported in more detail in gll|28|,|29|,|30|. 



3.2. The gravitational time delay 

The Lorentz-invariant relativistic time delay equation, generalizing the static Shapiro delay |3^ . can be obtained 
directly from Eq. H3.10|l . We note that the phase of electromagnetic wave, emitted at the point Xq = (ctQ^Xo) 
and received at the point — (ct,x), remains constant along the wave's path ^3. ■25„ .31.]. Indeed, if A is an affine 
parameter along the path, one has for the phase's derivative 
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which means that ip (a;" (A)) =const in accordance with our assertion. Equating two values of the phase at the points 
Xq and and separating time and space coordinates one obtains from H3.10|l 

1 , , 2GM, 1 - c-ifc -vj , ^ , , 

t — to ^ -k ■ (x — Xq) — hi [r — k ■ r) + const , (3-14) 

where 'const' denotes aU of the constant terms and the rclativistic terms taken at the emitting point, the distance 
r — x — xj{s), r = \x — xj{s)\, and the retarded time s is defined by the gravity cone Eq. (|2.7() . The Lorcntz-invariant 
expression for time delay (|3.14|l was derived first by Kopeikin and Schafer who solved equations for light geodesies 
in the gravitational field of moving bodies. Klioner |33l | also obtained this expression by making use of the Lorentz 
transformation of Jupiter's coordinate and the light vector in the Shapiro delay (which is implicitly equivalent to 
transforming the gravity and electromagnetic field equations) from a static to moving frame. Notice that in GR Eq. 
H3.14(l describes a hypersurface of both the light and gravity null cones along which electromagnetic field from the 
quasar and the gravitational field from Jupiter are propagating respectively (see Fig. P). The light cone intersects 
with the gravity null cone given by Eq. H2.7|l at the point a;" at which the propagating Jupiter's gravity field reaches 
the hypersurface of electromagnetic phase at time t when Jupiter is located at its retarded position xj{s) coinciding 
with the origin of the retarded gravity null cone. Additional details are shown in appendix where we used a bi-metric 
theory of gravity. In this theory the light null-cone does not coincide with the gravity null-cone as it was in GR, and 
it implies that light can not propagate on the hypersurface of the gravity null cone H2.7|) . Hence, the retarded position 
of Jupiter measured in Jovian experiment is due to the finite value of the fundamental speed of gravity and it sets 
the upper limit on the speed of propagation of gravitational waves which can propagate either on shell of the gravity 
null-cone H2.7|l in case of GR or inside it in the case of the theory of gravity with massive gravitons [s^ . 



3.3. The aberration of gravity 

Let us start with the light geodesies Eq. H3.9|l . Taking partial derivative from the metric tensor recasts Eq. H3.9() 
to the following form 



where |15||1 

= d^-TR = — - , (3.16) 

and we have neglected Jupiter's acceleration. 

The vector r" is null due to the gravity null cone Eq. (|2.t)l) . but the vector is space-like, because ricn" = +1. 
This occurs because the retardation of gravity effect in the first term in the right side of Eq. H3.1(j|l (present in 
the coordinate a;j(s) of Jupiter through the retarded time s) is compensated by the second term depending on the 
four-velocity of Jupiter which is a time-like vector, Mqu" = —1. This velocity-dependent term present in the 
gravitational force is associated in GR with the rclativistic effect known as aberration of gravity because it 
describes an aberration- like change in the direction of the gravitational force (that is, from the null- like vector r^/r^j 
to the space-like vector n") caused by motion of the gravitating body with respect to the reference frame used for 
data reduction. 

The compensation of the retardation of gravity by the aberration of gravity in the expression for vector can be 
demonstrated using a Taylor (post-Newtonian) expansion of a;j(s) around time t. This gives 

r" = i?" _ (s _ = i?" + , (3.17) 

ds c 7 

TR = -M^i?^ + — , (3.18) 

where terms depending on Jupiter's acceleration have been neglected, i?" = — x°'j[t) = (0,a; — xj(i?)) is a purely 
spatial vector lying on the hypersurface of constant time r = \x — a;j(s)|, and the retarded time Eq. H2.7|l has 
replaced s — t with r/c where c in this case is understood as the fundamental speed of gravity. Substitution of Eqs. 
(|3.17() . H3.18|l into Eq. H3.16|l . and reduction of similar terms, show that retarded term proportional to s — t cancels 
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with the gravity aberration term proportional to u" and Eqs. H3.16() and (|3.18|) are reduced to 



(3.19) 
(3.20) 



This result, shown elsewhere jlOj . reveals that (in the first approximation with respect to the four- velocity u") the 
gravitational acceleration of a light particle (photon) is directed toward the present position of Jupiter Xj(t) separated 
from the photon by a space- like vector R". At first glance, this might imply that the aberration of gravity cannot be 
measured in the linear vj/c approximation beyond the Shapiro delay. This point of view is advocated in [EIESBIS- 

However, it is crucial to understand that the VLBI experiment does not measure the gravitational acceleration (the 
right side of Eq. H3.15|l ) of radio wave (light particles) which is not VLBI-measurable quantity. What VLBI measures 
in reality is the electromagnetic phase of the radio wave from a quasar and direction of its propagation, which do 
depend on the retarded position of Jupiter xj{s) in contrast to the gravitational acceleration. In order to see how it 
happens let us integrate the inverse distance l/r/j along the light-ray trajectory parametrized by the affine parameter 
A. It results in 



dX 



1 ds 



1 



7 kar° 



ln(~fc„r") 



(3.21) 



which can be easily confirmed by direct differentiation of the gravity null cone Eq. (|2.6|l taken on the trajectory of 
the radio wave signal a;" = a;" (A). Using Eq. H3.21|l for the integration of the hght geodesic Eq. (|3.15|l gives 



2GMj k^uf" 



ka + (/CixU") — 

rR 



(3.22) 



which coincides with the expression for Ka ~ da'-p^ derived by differentiation of the electromagnetic phase 13.10|l . 
The subsequent post-Newtonian expansion (|3.17|) of r" in the function fear", that enters solution H3.22|l for and 
that H3.10|l for the phase 1^9(2;"), gives terms which depend on the aberration and the fundamental speed of gravity 
already in linear order ~ vj / c. This result is confirmed in j21„ 28] using general relativity parametrized by a single 
speed-of-gravity parameter Cg characterizing all non-stationary effects of gravitational field. 



3.4. Why can VLBI determine the retarded position of Jupiter? 

The gravitational force exerted by moving Jupiter on any mass (including radio photons) is a spatial vector directed 
in the first approximation toward the present position of Jupiter as shown in Eqs. (|3.15|) and (|3.19|) . A question arises: 
"How can the gravitational force be a spatial vector but the phase of electromagnetic waves, measured by VLBI, be 
gravitationally affected from null direction given by Jupiter's retarded position?" The answer to this question arises 
from the nature of gravity and radio waves which both move along null rays and their wave fronts form a null 
hypersurfacc (null cone). 

The gravitational 4-force on the right side of Eq. (|3.15|) has a Newtonian analog, and for this reason this force is a 
space-like vector directed to the present position of Jupiter x'^{t). Integration of the 4-force along particle's trajectory 
gives a vector of a linear momentum of the particle which is tangent to particle's world line. In the case of a slowly 
moving particle, its linear momentum is a time-like vector. However, radio photons move along light geodesies and 
their linear momentum is a null wave vector j^jE^. Hence, the integration of the light geodesies (|3.15|) along the 
photon trajectory must produce the null vector K'^ = fc" + 5k°' of the photon, where fc" is an unperturbed wave 
vector of the photon and (5fc" , given by Eq. (|3.22|) , is its gravitational perturbation caused by the space-like vector of 
the gravitational force in the right side of the equations of motion H3.15|l . 

The null vector K°' is a gradient of the radio phase with K°' = g"^dj3(p and the phase ip{x") is a null hypersurface. 
The gravitational field of Jupiter also propagate along null characteristics defined by vector = — x"{s), where 
Xj (s) is a retarded position of Jupiter taken at the retarded time s defined by the gravity null cone Eq. (|2.7|) . For this 
reason, the characteristics of the gravitational field of moving Jupiter form a null hypersurface, gravitational phase 
ipgrix"), as well. Radio photon interacts with the gravitational force at point x". Therefore, two null hypersurfaces 
iflx") and ipgrix") intersect at this point and form a null sub- hypersurface such that both vectors K"' and r" 
belong to it, and the gravitational perturbation of the radio phase, ip(x°'), is defined by the structure of this null 
sub-hypersurface. 

There are only three scalars that can be formed from the three vectors fc", r", and u" (Jupiter's velocity): kau", 
rR = —Uar°', and kar", and the (weak) perturbation of the radio wave phase can depend on a linear combination 
of some functions having these three scalars as arguments. The scalar r^j, however, can be eliminated because the 
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partial derivative of the phase 9q<p must be a null vector, whereas the partial derivative d'^rn = n" (see Eq. (I3.16() ') is 
a space-like vector which does not lie on the null hypersurface. The scalar k^u" is clearly a constant, if one neglects 
Jupiter's acceleration, so that d^{kaU°') = 0. Thus, term kau" has no influence on the propagation direction of 
the photons. We conclude that the only significant scalar combination defining the gravitational perturbation of the 
electromagnetic phase ip is simply k^r"', which depends on the retarded position Xj{s) of Jupiter taken at the retarded 
instant of time s defined by the gravity null cone Eq. (|2.7|l obtained as a retarded solution of the gravity field Eqs. 
(|2.3|) . The exact solution of the eikonal Eq. (|3.8|l confirms the above argument, and shows explicitly that the radio 
phase ip is sensitive to the retarded position of Jupiter xj(s), whereas the gravitational force, acting on photons as a 
space-like vector, is sensitive to the present position of Jupiter xj{t). 

In summary, the solution of the Einstein equations shows that gravity force propagates from a moving massive body; 
hence the retarded position of the body, if the speed of gravity propagation is not infinite, is a natural consequence 
of the theory. But, the interaction of gravity in the dynamics of slowly moving bodies {v <C c) will not produce a null 
vector since their four-velocities are not partial derivatives of a null hypersurface, and the retardation of grav ity effect 
is relegated to higher orders of w/c in the solutions of the equations of motion for these bodies [sE EEj^- On the 
other hand, the phase (or time delay) of the radio waves which is measured by VLBI, is sensitive to the retardation 
of gravity effect (retarded position of Jupiter) at order v/c because in GR gravity has the same speed as light which 
means that the fundamental speed in the Maxwell and Einstein equations is one and the same. This prediction of 
general relativity was the goal of experimental testing in Jovian experiment. 



3.5. Fundamental speed of gravity as a fitting parameter 

Our theoretical formulation of the problem of propagation of light in a time-dependent gravitational field of moving 
bodies predicts that the retarded coordinate of Jupiter Xj{s) in all terms depending on gravity field in Eqs. (|2.5|) ~ (|3.1U|) 
is shifted from its present position xj(t) due to the finite time the gravity takes to propagate from a moving massive 
body (Jupiter) to the field point x". To quantify this prediction experimentally we have developed a co-ordinated 
fiber-bundle parameterization of both the Einstein equations and light geod esies, based on a (single) speed-of-gravity 
parameter Cg. This Cg-parametrized approach is explained in our papers [2l|,|2l| ^- This Cg -parameterization decouples 
the null characteristics of electromagnetic and gravitational field already at the first order terms in v/c in contrast to 
the PPN formalism by Nordtvedt and Will |Q| which starts to distinguish the light and gravity null cones only in terms 
of the second order in v/c [S^. We assert that on each bundle of the parametric general- relativistic space-time we are 
allowed to substitute Eq. H2.7|l with the following description of the characteristics of the gravity field equations 

Sg ^ t ~ —\X - Xj{Sg)\ , (3.23) 

where Cg is a speed of gravity parameter running formally from Cg = oo to Cg = c. Again, we assert that in the papers 
|2lL l28j we are not using Cg ^ c in. the sense of an alternative theory of gravity. We work with various values of Cg in 
the fiber-bundle space-times of GR from which only one coincides with physical space-time when Cg = c. Parameter Cg 
must be considered rather as a label marking matrix of the Lorentz transformation of the gravitational field variables 
[28j | , making it clearly different from the matrix of the Lorentz transformation of electromagnetic field which depends 
on the speed of light c. For this reason Cg is not a measure of the PPN parameter ai which is associated with preferred 
frame- velocity tests of GR, described in ^2^] (see misconception 7) . 

Values of Cg less than c (although such values can be formally obtained from the fitting procedure) are strongly 
restricted by observations of cosmic rays. Indeed, if Cg were significantly less than c, the cosmic rays would radiate 
energy mostly in the form of gravitational Cherenkov's radiation, which contradicts the observations . Eq. (|3.23|) 
is valid up to the post-Newtonian terms of second order in the sense that it is obtained as a retarded solution of the 
Einstein equations where all time derivatives contain a single speed-of-gravity parameter Cg plL |2^ . The value of the 
parameter Cg in our model of the experiment was found by minimizing the phase residuals Sip = |<Pobs — ^csdix°')\, 
where (/?obs is the observed value of the phase, and 

(Pcal(x") = (^0 + fccx" + (fc„w") In (-fcap") , (3.24) 

which is the model function for the phase obtained from Eq. (|3.1(J|I after mathematically legitimate replacement of 
r" with = - x°^{sg) mil. 



Appendix of the present paper gives another theoretical parametrization of the speed of gravity versus the speed of Hght 
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We have introduced the speed-of-gravity parameter Cg to Eqs. (|3.23() and H3.24|l in such a way that eikonal Eq. 
(|3.8|) does not change the physical identity of the gravity null cone Eq. 1)2. 7|l which treats the retardation as due to 
the finite speed of the gravity field. The solution to the Einstein equations remains the same as in Eq. 1)2. 5|l with s 
replaced with Sg '2l', '28*1. This Cg speed-of-gravity parameterization of the Einstein equations demands that the speed 
of the gravitational field is the same in both the wave and near zones of an isolated system emitting gravitational 
waves. Our experimental goal was to measure Cg in the near zone, and to determine ii Cg — c predicted by general 
relativity. As we noted previously, the experiment can not measure the speed of gravitational waves in the wave 
zone directly but the measurement of Cg in the near zone definitely imposes the upper limit on the speed of their 
propagation everywhere. 

The retardation of gravity is unobservable in the solutions of the post-Newtonian equations of motion of self- 
gravitating bodies (a binary pulsar and the Earth-Sun, for example) in terms of order v/c and v^/c^ beyond the 
Newtonian law of gravity because of a perfect cancellation of the retarded and velocity-dependent terms present in 
the gravitational force between the bodies [3^|3^|33. However, the present paper (see also PjII^HoI) shows that 
the retarded terms are not canceled in the solution l|3.1()|l of the eikonal equation by the velocity-dependent terms. 
The phase ip has an argument /cqt" which does not depend on the velocity of Jupiter explicitly. The post-Newtonian 
expansion of fc^r" around the present time t generates velocity-dependent terms which are physically explained by 
the finite speed of gravity. 

Our Cg parametrization |^ assumes that the speed-of-gravity parameter Cg appears in the equation for the 
time delay t — Iq taken by light to propagate from Xq to x in the following form 

where /c" is a wave vector of light, v is its frequency, and fe — (fc') is a unit vector along the direction of its unperturbed 
propagation. Integration of Eq. (|3.25() leads exactly to formula H3.14|l with the retarded time s replaced with Sg defined 
in Eq. (|3.23|l of this section. It is clear that in case Cg = oo the Jovian deflection experiment would not measure 
the integral from the time derivative of the metric tensor shown as the second (gravito-magnetic) term in the right 
side of Eq. (|3.25|) . However, we have measured Cg ~ c with the precision 20% and, thus, proved that the experiment 
is sensitive to the gravito-magnetic dragging of light generated by the partial time derivative of the metric tensor 

MM 



,i dh^^{T,x) 1 dhf,i,{T,x) 



dr 



dT 



(3.25) 



3.6. The determination of the speed of gravity 



The speed-of-gravity parameter was determined from the experiment in the following way. Let us introduce two 
angles 9 = 8(i) and 6 = 0{sg). The angle O is between the unit vector k characterizing direction of propagation 
of the light ray and the unit vector p = R/\R\ while the angle 6 is between vector fe and I = l{sg) — p/p. Here 
R = X — xj{t) connects the present position of Jupiter xj{t) at the time of observation t and the point of observation 
X, whereas p = x — xj{sg) connects the retarded position of Jupiter at the retarded time Sg — t — p/cg and the point 
x (see Eq. H3.23|l ). By definition cos 8 = fc • p, and cos 9 — k ■ I (see figure 

From the definitions of vector l{sg) and the angle 0{sg), we can obtain their post-Newtonian expansions in the near 
zone of the solar system by using Eq. H3.17|l . Notice that l{t) — hms^^f i(sg) = p, and 9{t) — liuis^^t 9{sg) = 6. 
The argument of the logarithm in Eq. H3.24|l for (pcai is 



-k^p- = ^ 
c 



1- cos 6l(sg) . (3.26) 
The post-Newtonian expansion of the model function of the electromagnetic phase yields 



<Pcai ^^0 + kax'' + {kc.u") ln|l - k ■ [p + -p X {v J xp)]} + f 4) : (3-27) 

C I. Cg ^ J 

where the constant and Inp terms have been neglected because they are smaller than the other relativistic terms, and 
are much smaller than the experimental angular resolution. 

The term in square brackets of the logarithmic function describes a small displacement of vector p to the retarded 
direction 

l=p+—px(vjxp). (3.28) 
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This equation is independent of the electromagnetic wave vector k and describes a change from the present (ephemeris) 
position of Jupiter given by vector p to its retarded location given by vector I (see figure This displacement is 
caused by the propagation of gravity force in the near zone of Jupiter with finite speed, and is associated with the 
near-zone aberration of Jupiter's gravity field. 

The determination of ipo, the unperturbed direction to the quasar fc, and coordinates of Jupiter in ipcai can be 
obtained to sufficient accuracy using the JPL ephemeris of the solar system objects with respect to the quasi-inertial 
reference frame, defined by the International Celestial Reference Frame grid of quasars in the sky This allows the 
measurement of the retardation of gravity incorporated in Eq. (|2.7|) for retarded time s. 

As we emphasized in previous sections, the aberration of gravity H3.28|l is unobservable in the case of the gravitational 
interaction between two slowly- moving bodies |lCll| . However, the gravitational perturbation of the electromagnetic 
phase (not the acceleration) of the radio wave was measured in the experiment in order to extract information about 
the direction to the position of Jupiter which, in addition to the well-known (radial) Einstein light deflection, has 
gravitational drag of the radio wave (tangential deflection) by means of the gravito-magnetic force produced by its 
orbital motion. This gravito-magnetic dragging is equivalent to the aberration of the gravity force. 

At the time of the closest approach of Jupiter to the quasar on 2002 September 8, the angle O ^ 3.7' = 10~^ rad, 
so that 1 — k ■ p = (1/2)9^ -t- 0(9^). On the other hand the post-Newtonian parameter vj/c = 4.5 x 10~^ <C 6. 
Hence vj/{cQ) <C 1, and the aberration of gravity post-Newtonian Eq. (|3.27|l can be expanded with respect to a 
small parameter vj/{cQ). It yields equation 



AGMj 



ine 



k-vj -{k- p){p ■ vj) 



0(4 

^9 



(3.29) 



which was used to estimate the magnitude of the aberration of gravity effect ,2j described by the second term in square 
brackets. The magnitude of the aberration of gravity is inversely proportional to the speed-of-gravity parameter Cg 



which has been measured in this experiment as Cg = 1.06 ± 0.21 in geometrized units Higher angular accuracy 
is needed in order to detect TT gravitational waves and measure their speed directly in the wave zone with existing 
pof or future 0| technologies of the gravitational- wave detectors and/or other technique [H ]. 



4. COMMON MISCONCEPTIONS ASSOCIATED WITH THIS EXPERIMENT 
4.1. Misconception 1: The gravitational field of a uniformly moving body is static 

This misconception is associated with the following reasoning. The gravitational field of a uniformly moving 
body can be obtained by solving Einstein's equations in the static frame of the body with its subsequent Lorcntz 
transformation to the moving frame. Since the gravitational field in the body's reference frame is static and not 
propagating, this "non-propagation" property of the gravitational field will be preserved after making the Lorentz 
transformation to the moving frame (see, for example, H^])- 

Any field is called propagating if it is obtained as a solution of the wave-type equation. A field is static if, and 
only if, it is obtained as a solution of elliptic-type equation. But the linearized Einstein's equations in the moving 
frame are not of the elliptic-type, they contain time derivatives and are essentially of the wave Eqs. (|2.3|l . The causal 
(retarded) solution of these equations describes the gravitational field of a uniformly moving body which depends on 
time. This dependence is represented in the form of a causal gravity wave, Eq. H2.5|) . propagating from the moving 
body (Jupiter) to the field point x" along the hypersurface of gravity null cone (|2.7|) . The gravity null cone structure 
is hidden in the static solution of Einstein's equations in the rest frame of the body. But, a Lorentz transformation 
from the static frame of the body to a moving frame preserves the structure of the gravity null cone, as Einstein's 
equations insist, but now makes the cone characteristics measurable by observing the retarded position from which 
the gravity field of the moving body deflects light (see Figs. 1-3). As discussed previously, the null-cone gravity-field 
characteristics become measurable in the moving frame because the time derivatives of the Einstein equation are 
normalized to the fundamental speed (of gravity), and these terms produce the retarded gravitational field, which is 
not static. 



* Notice that this expansion is vah d un til the relationship v/c < Q is hold. It is incompetent to take the limit — > for fixed value of 
d/c as mistakenly done by Asada l39l 
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4.2. Misconception 2: The propagation of the gravitational field of a uniformly moving body can not be 

observed 

The gravitational force acting on a test particle moving in the field of a uniformly-moving massive body is directed 
in the first approximation to the present position of the body in accordance with Eqs. H3.19|l and (|3.2U|) . Thus, it 
is assumed that the propagating property (retardation) of the gravity field of a uniformly moving body cannot be 
observed (i^. 

This is true if the trajectory of the test particle is the only measurable quantity. The normal use of relativistic 
celestial mechanics generally deals with such problems as the motion of binary systems and solar system objects 
which velocities are much smaller than the fundamental speed of gravity. However, this is not true if the phase of an 
electromagnetic wave is measured as its front propagates in the field of a uniformly moving body. As GR predicts and 
this is shown in section I? 3. 41 the phase depends on the retarded position of the body. This retarded position is caused 
by the finite value of the fundamental speed of gravity in the gravity null cone Eq. (|2.7|l ; hence the fundamental speed 
for gravity can be measured (see appendix for further details). 

4.3. Misconception 3: The parameter c of the Lorentz transformation is physically the speed of Ught 

The parameter c of the Lorentz transformation is generally associated with the speed of light [El^ilS' However, 
verbatim interpretation of c as a physical speed of light narrows the true meaning of the Lorentz transformation as 
describing the transformation property of Maxwell's equations only. Einstein's equations are also wave equations where 
the wave operator depends on the fundamental constant c which must be interpreted as the speed of propagation of 
gravitational field. Consequently, the parameter c in the Lorentz transformation associated with the Einstein gravity 
field equations must be related to the speed of gravity. Generally speaking, parameter c of the Lorentz transformation 
is just a fundamental constant of the Minkowski space which characterizes an ultimate speed of propagation of any 
massless field in vacuum. Hence, the name of the "speed of light" for the parameter c in the Lorentz transformation 
should not be taken literally when one interprets gravitational experiments dealing with non-stationary (propagating) 
gravitational fields. 

4.4. Misconception 4: The speed of gravity can only be measured through observation of gravitational waves 

The prevailing view is that the speed of gravity can only be associated with the propagation of plane gravitational 
waves decoupled from the system emitting these waves [43 . Such waves are generated if the second time derivative 
of the quadrupole moment of the system is not constant 3,"25"|; that is, the bodies comprising the system accelerate. 
These gravitational waves, which are type N according to Petrov's classification of gravitational fields JJe ^JIj 
are indeed generated if, and only if, the bodies accelerate. However, the gravitational field of a system emitting 
gravitational waves also has contributions from other Petrov's type fields - I, II, and III, which have in general 
relativity the same structure of the null characteristics as the N-type waves, but decay faster with distance from the 
system so that these near-zone fields can not be observed as freely propagating waves of type N 0|. Hence, any 
measurement of the null characteristics of the near-zone gravitational field provides us with the upper limit on the 
speed of propagation of gravitational waves of type N. 

Propagation of light through the multipolar gravitational field of an isolated astronomical system, emitting gravita- 
tional waves, in the near, intermediate, and wave zones of a system, have been already studied p^l45ll46i |. We showed 
that parameters of the light signal (wave vector, polarization, etc.) depends in general on the retarded value of the 
multipole moments of the system, as well as their first (velocity associated) and second (acceleration associated) time 
derivatives. The second time derivatives of the multipolar field dominate in the wave zone but plays no significant role 
in the near and intermediate zones where the non-radiative gravitational fields associated with the retarded value of 
the multipole moment of the system and its first derivatives prevail. This fields and their propagators can be studied 
in the light-ray deflection experiments when light moves in the near zone of the isolated system. Thus, a propagating 
gravitational field should not be associated only with the accelerated motion of the bodies when measurements are 
made in the near and/or intermediate zone of the system emitting gravitational waves. 

In summary, the propagation of gravity is a general concept which must be associated with the wave nature of 
gravitational field, and is not limited exclusively to plane gravitational waves of Petrov's N type. Thus, a uniformly 
moving body generates a retarded gravitational field which is associated with velocity of the body. This retarded 
velocity-field can be also obtained by making Lorentz transformation of the static gravitational field. This transfor- 
mation induces the aberration of the null characteristics of the gravity field which can be observed using the phase 
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of the radio waves and used to measure their speed without direct detection of gravitational waves, as explained in 
previous sections of this paper. 

The velocity- field is well-known in electrodynamics 0, l25l \4l\ when electromagnetic field of moving charge is 
calculated in terms of the retarded Lienard-Wiechert electromagnetic potentials. This velocity- field reflects the prop- 
agation of the electromagnetic field of a uniformly moving charge but it can be also calculated by making use of the 
Lorentz transformation. The same is true for gravity. 

4.5. Misconception 5: The retarded position of Jupiter is taken at the time of the light-ray closest approach 

The consensus model used to evaluate the deflection of light by Jupiter (egs. calculates the time t* of closest 
approach of the unperturbed light ray to Jupiter, determines the 'retarded' position of Jupiter xj{t*) at this time, and 
then calculates the deflection magnitude and direction. This interpretation of the retarded position of the light-ray 
deflecting body (Jupiter) sounds very plausible but the mathematical solution of the Einstein and eikonal equations 
reveals that the retarded position of Jupiter must be actually taken at the point Xj{s) on Jupiter's world line at the 
retarded time s defined by the solution of the gravity null cone Eq. I|2.7f) . 

General relativity unequivocally predicts that the deflection of the radio photon occurs the most strongly when 
Jupiter is at its retarded position Xj{s) on the gravity null cone (|2.7|l . In the case of a small impact parameter of the 
light ray to Jupiter, the time of the closest approach t* of the light ray to Jupiter is nearly equal to the retarded time s, 
Eq. 1)2. 7|) . caused by the finite speed of gravity. Thus, the retarded position xj{s) of Jupiter can be expanded around 
the time of the closest approach, xj{s)) — xj{t*) + {s — t*)vj{t*) + ... so that position of Jupiter xj{t*) taken at the 
time of the closest approach can be used in calculation of the light-ray deflection with a very good approximation. 
One should not however confuse the cause and effect - the time t* of the closest approach of the photon to Jupiter is 
not inherent to the theory at all, but comes from an approximation of small impact parameter of light-ray trajectory 
to the light-ray deflecting body. It does not mean that the retardation of the gravitational perturbation of light is 
caused by light propagation. Full treatment of the small and large impact parameter approximations in the case of 
gravitational Icnsing is given in [T^ ITM Isil . 

4.6. Misconception 6: Null characteristics of gravity field can be observed only at {v/cf' terms beyond the 

Shapiro delay 

This misconception arises from the post-Newtonian expansion of the metric tensor perturbation h^p in Eq. H2.5() 
made prior to solving the eikonal Eq. (|3.8|l . Indeed, the post- Newtonian expansion of /iq/3, using the Taylor expansion 
of the retarded variables depending on the retarded time s in Eq. H2.7|l around the present time t, leads to 

hap{s,x) = hap{t,x) + terms of order (v/c)^ , (4.30) 

where hap{t,x) are solutions of the elliptic- type equations (Newtonian potentials). Eq. H4.30() suggests that the 
retardation of gravity effect is postponed to the terms of order (v/c)^ because no terms of order v/c is explicitly 
present in this expansion 

The physical reason for the absence of terms of order v/c va Eq. H4.30|l is that the Lienard-Wiechert potentials 
depend on distance = —Uar" which incorporates both the retardation and aberration of gravity effects that are 
equivalent in terms of order v/c and cancel each other when the post-Newtonian expansion of the metric tensor is 
done. Light propagates in GR with the same speed as gravity. Lorentz-invariant solution of the eikonal equation 
incorporates the aberration of light which is equal to the aberration of gravity terms and mutually cancel each other, 
hence, leaving the retardation of gravity terms in the radio wave phase unchanged. That is why the post-Newtonian 
expansion of the retarded gravitational perturbation of the electromagnetic phase reproduces the aberration of gravity 
terms which we had observed in the Jovian deflection experiment via the retarded position of Jupiter 0] . 

Although hafj{t,x) does not contain terms of order v/c, the measurement of these Newtonian- like gravitational 
potentials by VLBI is indeterminate because the gauge invariance of the Einstein equations make the potentials 
coordinate-system dependent. There are gauges, like the ADM gauge 0, where 500 and goi components of the 
metric tensor become the solution of the elliptic-type equations which suggest that these metric tensor components 
do not propagate at all. But, of course, any gravitational influence of these potentials on measurable quantities is 
not instantaneous but propagates with the finite speed. This resembles the situation in electrodynamics where the 
Coulomb gauge leads to the elliptic-type equation for scalar potential but it does not mean that electric and magnetic 
fields propagate instantaneously [4^ . Hence, the potentials are not directly measurable quantities in the light-ray 
deflection experiments because they depend on the specific choice of gauge conditions which is a matter of convenience. 
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Only those observable quantities which are gauge-independent (like the electromagnetic phase) can be used to probe 
physical properties of the gravitational fields. 

Our method of analysis, given in this paper, consists of two part: first, to solve the problem of the propagation 
of light in the time-dependent gravitational field by making use of the Licnard-Wiechcrt potentials; and then to 
determine the post-Newtonian expansion of the solution in order to produce the physical interpretation in terms of 
different order in vjc. However, we will show below that if the post-Newtonian potentials (|4.3Q(I are used to solve 
the eikonal Eq. H3.8|) . the solution will contain terms of order vjc which can be represented in several equivalent 
forms — leading to difficulty in the interpretation of the solution. This difficulty can be resolved if the solution of the 
eikonal equation at higher orders of vjc is also determined in order to match it smoothly with the linear terms of 
order vjc present in the phase. 

For this experiment, the observable quantity is electromagnetic phase Lp which is a gauge-independent quantity 0. 
We can measure the impact of the gravitational potentials on the electromagnetic phase in the following manner. The 
post-Newtonian expansion of the electromagnetic phase is 

^(a;") = (^o + fcaa;" + ^^/^(a;") , (4.31) 

where (^o is constant and fc" is the wave vector of electromagnetic phase. Substitution of Eq. (|4.31|l to the eikonal 
equation and making use of expansion H4.3()|l yields the linearized equation for the phase perturbation i/; 

where Uj{t,x) — GMj/R, R = \x — xj{t)\, and fc' is the unit vector in the direction of propagation of light rays 
which are orthogonal to the hypersurface of the phase (p. Eq. (|4.32|) is a partial differential equation of the first order 
and its solution is 



V'(x") = 2GMj(^l - ^--^^ In +o(^^^ , (4.33) 



where p — k — c ^vj, and vj = dxj{t)/dt is Jupiter's velocity. 

The argument of the logarithmic function in Eq. I|4.33|) can be represented in two different forms. The first form 
is given bv [6Ll2^ 

pR-p-R^ [l--k-v)[R-K-R] , (4.34) 



c 



where 



The second form is 



K=-^k^-(kx(vjxk)]. (4.35) 

B C V / 



pR-p R = r-k r = -k^r"' , (4.36) 

where fc" is a wave vector of the radio wave, r" — {r,r), r = x — xj{s), r ~ \r{s)\, and the retarded time s is given 
by the gravity null cone Eq. H2.7|l . 

Eq. (|4.35l) looks similar in form to the aberration of light formula transforming the direction of propagation of 
the light ray k from the moving frame of an observer to the direction /C in the static frame of Jupiter. This is 
why Will 01 and Carlip Q interpreted the effect we have observed as the aberration of light (see next section and 
appendix). However, Eq. (|4.36() clearly shows that this aberration term should be interpreted as an aberration of 
gravity caused by its propagation because Jupiter's coordinate depend on the retarded time s defined through the 
solution of the Einstein (wave) Eq. H2.3(l . Any aberration, of course, is of order v/c but only one type of aberration of 
either electromagnetic or gravitational field can be unambiguously extended to terms of second order in v/c to keep 
Lorentz-invariance of the equations. Our calculations show that aberration term in Eq. (|4.35|l matches smoothly with 
the retardation of gravity effect H4.36|l in any post-Newtonian order while the aberration of light does not (see also 
|33| for more detail on the matching procedure). 
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4.7. Misconception 7: The Jovian experiment and the PPN parameter qi. 

The Jovian deflection experiment measures the aberration of gravity and its fundamental speed hmit. It tests for 
the first time the non-stationary gravito-magnetic property of the gravitational field, of order (v/c)^, generated by 
the time derivative of the metric tensor (gravito-magnetic displacement current). The experiment can also be viewed 
as a test of the Lorentz-invariance of the gravitational field and confirmation of its relativistic causal nature. 

Will 20] introduced the "preferred frame" PPN parameter ai which characterizes the violation of the Lorentz- 
invariance of the PPN metric tensor. However, this parameter presents only in goi component of the metric tensor 
if it is written down in the preferred frame. Hence, parameter ai characterizes, in fact, a strength of the physical 
coupling between the goi component of the metric tensor and the matter current given by Toi component of the 
stress-energy tensor. Subsequent Lorentz transformation to a frame moving with respect to the preferred frame with 
velocity w introduces a product aiw to goo and goi components of the metric tensor. 

Under the assumption that the cosmic microwave background radiation (CMBR) defines the preferred frame, several 
accurate experimental tests give ai < 4 x 10~^, which appears more accurate than our limit on (cg/c) — 1 < 0.27, 
both of which are measures of the experimental departure from GR. However, direct comparison between these two 
methods for testing Lorentz-invariance is not valid. First, the PPN formulation states that ai is not coupled with 
time derivatives of the metric tensor whereas the deflection experiment is a direct test of the magnitude of the time 
derivative of the metric tensor (see Eq. 13.25|l and 0,|2^). In GR the parameter ai = 0. Therefore, search for ai is 
an attempt to flnd an effect, which probably does not exist at all because GR is a valid theory, while our measurement 
of the fundamental speed of gravity discovers a real (and newly measured) physical effect - the aberration of gravity, 
which is a fundamental property of gravity predicting that electromagnetic signals are deflected by gravitational fleld 
of a moving gravitational body from its retarded position with respect to observer. 

Secondly, the current limit on the parameter ai strongly depends on our guess about the velocity w of the solar 
system relative to a "preferred frame" which was assumed to be the CMBR. But, we do not know if it defines the 
preferred coordinate system in the sense forbidden by special relativity; it is just a convenient frame to describe the 
global isotropy of the cosmological black-body radiation. Furthermore, future observations of the relic gravitational 
wave background (GWB), for example, may lead to the frame somehow moving with respect to the CMBR frame. 
Since the GWB was formed in the very early universe, long before the CMBR decoupled from matter, it would be 
more likely to associate it with the preferred frame. However, neither the global topological isotropy of CMBR nor 
GWB may be related to the local Lorentz-invariant isotropy of the space-time itself. Finally, modern multi-connected 
cosmologies challenge the cosmological Copernican Principle, and indicate that various possibilities for the global 
preferred frame may not coincide with the CMBR and GWB frames Thus, the limit on ai, established by |23| 
depends on numerous assumptions about properties of a preferred reference frame and other cosmological principles 
which are to some extent ambiguous concepts. It seems to us that only GP-B experiment will be able to set up a 
limit on ai parameter which will not be sensitive to the choice of the preferred frame |50j |. 

5. OTHER INTERPRETATIONS 

Four alternative interpretations, which disagree with our interpretation and among themselves, of the Jovian deflec- 
tion experiment have been proposed (see section^. We have analyzed fully the alternative interpretations elsewhere 
HE HOl I but summarize the main arguments here. 

The main reason for the disagreements between the four other interpretations is their use of the truncated linearized 
expansion of the Einstein equations which is not capable to study the Lorentz-invariance of the gravity fleld force in 
the light-ray deflection experiments in higher orders in v/c. The truncated linearized theory is not sufflcient to give 
the proper geometrical picture of the observed quantity in the time delay of light which is, in fact, a Minkowskian 
dot product $ — —kar" between two null vectors fc" and r" (see flgure^ entering the gravitational perturbation of 
the electromagnetic phase (p in Eq. H^.lOfl . When Einstein's theory of the gravitational time delay is truncated at 
linear v/c order terms, the simple geometric picture of the experiment is easily missed and the interpretation of the 
experiment becomes obscured. 

Win H (see also 0) introduced a set of alternative theories of gravity with a propagation speed Cg which could 
differ from the speed of light c. However, the use of the Post-Newtonian Parameterized (PPN) approximation, which 
is wide-spread, does not keep the geometrical relationship between the Christoffel symbols (light geodesic equations) 
and Ricci tensor for different values of Cg. Moreover, the PPN formalism transforms the Christoffel symbols from one 
inertial frame to another with the matrix of the Lorentz transformation for electromagnetic fleld p3 | . This makes no 
sense since the Christoffel symbols are deflned purely in terms of gravitational potentials which must be transformed 
in accordance with the group of the Lorentz transformation of gravitational fleld depending on the parameter Cg. In 
addition, the PPN gravity fleld equations with Cg are not derivable from a Lagrangian of any known theory of gravity. 
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Hence, they do not provide us with a complete physical description of the gravitational field in case when the speed 
of gravity Cg ^ c making the PPN approach unreliable for analysis of experiments in time-dependent gravitational 
fields. This may explain why Will's use of the PPN approximation to analyze relativistic corrections to the Shapiro 
time delay confuses the aberration of light and the aberration of gravity (see appendix for more detail) . 

As discussed in the previous section, Will's formulation does predict an 'aberration' deflection of order tj/c, Eq. 
(|4.35(l , although Will attributes this aberration to propagation of "light" . The difference between his interpretation 
with ours can be resolved by finding explicit terms of the second order in vj/c in the solution of eikonal Eq. (|3.8|l . 
and then matching them smoothly with the solution (|4.33|) . This matching shows [2ll | that the correct interpretation 
of the retarded position of Jupiter in the gravitational perturbation of electromagnetic phase in Eq. H4.33|l must be 
associated with the null characteristics of gravity field. In other words, had the exact aberration-of-light formula been 
used to account for the second order terms in Eq. (|4.35|) . it would not produce all of the second and higher order 
terms in the solution of the eikonal Eq. 13.8|l . whereas the retarded time Eq. (|4.3()|l gives all terms of second and 
higher orders in the Lorentz covariant solution of the eikonal Eq. H3.8() as shown by the calculations given in the 
present paper and in |2l,ll0)ll3- Thus, the correct meaning of the linearized Eq. (|4.35(l should not be associated 
with the aberration of light but with the retardation (aberration) of gravity caused by the orbital motion of Jupiter 
[2^ and associated with the finite value of the fundamental speed of gravity ^ . 

Asada Q suggested that the origin of the retarded time Eq. (|2.7|l is associated with the propagation of the radio 
waves from the quasar in the field of static Jupiter. However, this equation originates from the Lienard-Wiechert 
solution of the Einstein equations (|2.3|) , which describes the propagation of gravitational field and is not associated 
in any way with the propagation of radio waves from the quasar. Propagation of radio waves is described by the 
time delay Eq. (|3.14|l which describes the light cone, while the retarded time Eq. (|2.7|l defines characteristics of the 
gravity field null cone (see figure Calculations displayed in the appendix of this paper obviously demonstrate that 
Asada's interpretation 5] is illicit since it confuses the null characteristics of light and gravity. 

A fundamental flaw in Samuel's Q interpretation was his assumption that the direction to Jupiter was directly 
measured by VLBI network in the deflection experiment so he also confused the null directions associated with 
characteristics of gravity and radio waves |29|. Even at the minimal separation of 3.7' of Jupiter and the quasar 
on 2002 September 8, no radio emission from Jupiter was obtained with the VLBI observations. The null direction 
to retarded position of Jupiter, depending on the value of the fundamental speed of gravity, was measured from 
the gravitational time delay Eq. (|3.14|) and compared with that known extremely accurately from the JPL solar 
system ephemeris depending on the speed of light (see appendix and Fig. Therefore, minimization of the residual 
electromagnetic phase with the speed of gravity taken as a fitting parameter allowed us to measure it with respect to 
the speed of light. 

Pascual-Sanchez 9j stated that the retarded position of Jupiter measured in the experiment is the R0mer delay 
of light coming from Jupiter to observer and that this R0mer delay has nothing to do with the gravitational time 
delay of light. Eq. 1)2. 7|) for retarded time s looks indeed similar to the R0mer delay but this delay is for propagation 
of the gravitational field from moving Jupiter to the observer at Earth as independently proved in the appendix of 
the present paper and in It does not describe the R0mer delay of light but gravity which takes finite time to 
propagate from moving Jupiter to observer. Had the R0mer delay of gravity been different from that for light, Jupiter 
would deflect light not from its retarded but from some other position on its orbit. The experiment disproved this 
assumption confirming that the R0mer delay of gravity and light are equal with accuracy of 20% by comparing two 
positions of Jupiter on its orbit obtained independently from the gravitational deflection of quasar's light and JPL 
ephemeris. 

6. SUMMARY 

We have developed a Lorentz-invariant general-relativistic theory of light-deflection experiments associated with 
moving bodies in the solar system. The formulation is based on the retarded Lienard-Wiechert solution of the Einstein 
equations, and explicitly calculates the effect on the phase of the radio wave from the quasar caused by the finite speed 
of gravity which is a primary fundamental constant in Einstein's field equations. Our post-Newtonian expansion of the 
retarded coordinate of the light-ray deflecting body shows that the retardation of gravity effect after making its post- 
Newtonian expansion, is present already in linear aberration-like terms of order v/c beyond the Shapiro delay; hence, 
we interpret the v/c term as caused by the aberration of gravity in the expression for the gravitational force (|3.15|) 



AS we mentioned the effect we have observed can be also explained in terms of the gravito-magnetic dragging of the quasar light ray 

mm 
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exerted by Jupiter on the photons. We have used the concept of the gravito-magnetic field caused by translational 
motion of Jupiter as an aid in understanding the gravitational phenomena associated with this experiment. 

The VLBI measurements confirmed that the deflection is associated with the retarded position of Jupiter, with 
the speed of gravity equal to that of light to 20% accuracy and this retardation is not a preferred frame effect 
associated with the PPN parameter ai. Other contradictory interpretations of the experimental results are associated 
with: (1) an analysis of the Einstein equations which are not fully Lorcntz- invariant; (2) misconceptions about the 
experimental method and analysis; (3) misunderstanding of the nature of the interaction of the gravitational field 
from moving objects with the phase of the quasar radio signal which is a null hypersurface. According to GR the 
phase is influenced in all orders of v/c by the propagating gravitational field of Jupiter and this perturbation reveals 
itself in the linearized order of v/c as the aberration of gravity force of Jupiter due to the finite speed of gravity. 

We thank the Eppley Foundation for Research (award 002672) for support. 



APPENDIX A: PROPAGATION OF LIGHT IN A BI-METRIC THEORY 

In this section we assume that the Einstein theory of general relativity is valid with the time coordinate x'^ — ct 
where c denotes the fundamental speed of gravity. Let us assume that light propagates in vacuum with the speed 
ci less than the speed of gravity c. Effectively, this means that vacuum in the presence of the gravity field can be 
considered as a fully transparent medium which has a constant refractive index e > 1. Propagation of light in such 
vacuum can be treated on the basis of the Maxwell electrodynamics in medium. The geometric optics limit of these 
equations has been elaborated by Synge |5l| and we employ his approach. 

According to Synge |51j light propagates in the medium along light geodesies of an optical metric defined by 

gaf3 = gap ^ ' 

where gap is the gravitational metric given by Eqs. H2.1|l . H2.5() . and is a vector field associated with four- velocity 
of the medium. In our case it means that if the speed of gravity and light are not equal there must exist a preferred 
frame in which the speed of light is isotropic while in all other frames it is anisotropic. Notice that in the preferred 
frame = (1, 0, 0, 0) and the speed of fight c; — c/e < c. 

There are other bi- metric theories of gravity (see, for example, 0ill3|) with two metrics predicting different speeds 
of propagation for light and gravity. In such theories the four- velocity V"' is identified with a long-range vector field 
spontaneously violating the Lorentz-invariance of gravity. Those bi-metric theories lead to the geodesic equations of 
light also propagating along the null cone of the optical metric . 

Let us assume for simplicity that both the refractive index e and the vector field V" are constant ^ . This model 
is simple but self-consistent and sufficient to apparently distinguish the speed of light ci from the speed of gravity c. 
Equations of our bi-metric model and their solution are given in the rest of this appendix. 

Light rays in Synge's approach [sif are defined by a covariant equation for electromagnetic phase (eikonal) which 



reads 



r^d^^d,^ = , (A.3) 

which generalizes Eq. H3.8|l for propagation of light in the case of c; 7^ c. Assuming that unperturbed solution of Eq. 
(jA.3|l is a plane wave we can write a general solution of this equations as follows 

(/j(2:") = </Jo + fcaa;"+V(a;") , (A.4) 

where ka is an unperturbed (constant) wave co- vector of the electromagnetic wave, and 'ip{x) is a relativistic pertur- 
bation of the eikonal generated by the metric tensor perturbation hap defined in Eq. I|2.5|l . Substitution of Eq. I|A.4|) 
to HA.3|I yields 

fj^^kak^ = 0, (A.5) 
V"^ka-^ = ^h'^^kakp , (A.6) 



More general case can be solved as well but it significantly complicates solution of equations and will be omitted 
Lagrangian-based approach of an alternative bi-metric theory of gravity Q leads to the same equation. 
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where the unperturbed part of the optical metric ljA.l|l . (|A.2p is defined by 



- (e^ - 1) VV^ , (A.7) 

j7a/3 = Val3 + (^1 " ^) Vc^Vp . (A.8) 

Let us introduce a Hght-ray vector a" defining direction of propagation of the unperturbed light ray. One has [sij 

cr" = f^^l^kf} = fc" - (e^ - 1) (V^^kp) V" , (A.9) 

which implies that 

fjapa'^a'^ = , and kacr" = . (A. 10) 

Vector cr" points to the direction of propagation of a light ray from a radio source (quasar) to observer (see Fig. 
Making use of vector a" simplifies Eq. HA.6|I and reduces it to the following form 

a'^pL ^ i/i„.a"(T^ . (A.ll) 

The unperturbed characteristics of the eikonal Eq. (|A.11|I are straight lines (light rays) parametrized by the affine 
parameter A in such a way that 

:^=a"^. (A.12) 

Integration of Eq. HA.12|I by making use of the unperturbed characteristics is straightforward (see, for example, (52i) 
and is given as follows 

. _ 2GM + (1/2) _ 

where we have defined 

= cr^; +(7xu" , (A. 14) 

al = P'^y , (A. 15) 

a± - (a^La'Tj)'/' - {Pa.pa'^aP)'^' , (A.16) 

and 

Pal3 = ??Q/3 + UaUp , (A. 17) 

is the operator of projection on the hyperplane orthogonal to the four-velocity of the light-ray deflecting body so 
that it obeys to equation PaijP^ — Pa-y- It is easy to confirm that Eq. (jA.13|l is solution of Eq. (jA.ll|) by observing 
that 

-il„(-V-..)^-fi, (A.18) 

where rjj = —u^r'^, and Eqs. H3.11|l . (|3.12|) have been used. Eq. (|A.14|) allows us to recast the argument of the 
logarithm in Eq. (jA.lSjl as 

^ar" = fTjr„ - d^rfl . (A.19) 
It is worth noticing that both vectors and r" are null vectors of the unperturbed gravity metric rjap, that is, 

Vcpl'^l" = 0, (A.20) 
Tlcpr'^rP = 0, (A.21) 

which are consequences of definitions given by Eqs. ljA.14|l and H2.6() . Thus, neither nor r" belong to the null cone 
of the optical metric gap (see Fig. (S)). This point has not been understood by researchers whose interpretations of the 
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Jovian deflection experiment differ between themselves. In particular, Asada |5'| insisted that the retarded position of 
Jupiter is on the null cone of the optical metric. Our calculations show that this Asada's assertion is erroneous. 

Solution (|A.4|I . I)A.13(I for the electromagnetic eikonal in the bi- metric theory should be compared with a similar 
solution H3.10|l for the case of propagation of light in general relativity where gravity and light have the same null 
cone. The reader can see that the null directions l|2.7|l of the gravity metric g^p enters explicitly the gravitationally 
perturbed part of the eikonal (|A.13|I in the bi- metric theory in the form of a dot product /qt" in the argument of 
the logarithm, where r" is the null distance of the metric gaf} between observer and the light-ray deflecting body 
®. This term replaces the argument kar°' of the logarithmic term in the eikonal (|3.1l)|l in general relativity where 
ci = c. A remarkable fact is that both P and r" are null vectors of the metric gaf3 describing propagation of gravity. 
Consequently, gravitational light-ray deflection experiments in the field of moving bodies are sensitive to, and can 
measure, the null characteristics of the gravity metric by observing propagation of light that moves along the null 
characteristics of the optical metric gap- This fully agrees with our general relativistic interpretation of the Jovian 
experiment 

Analysis given in this appendix provides a firm support to our statement 0, Q that the retarded coordinate of 
Jupiter Xj(s) measured in the Jovian experiment is associated with the finite speed of gravity but not light. The 
statement that the Jovian experiment measures the light-cone effects associated with the optical metric gap as claimed 
in Hi, 0,11 Hm is not valid. 



The amplitude of the eikonal also changes leading to the deviation of the PPN parameter 7 from its general relativistic value 7 = 1. 
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FIG. 1: The two null cones related to the experiment are shown. The gravity cone is a solution of the Einstein equations and 
describes null characteristics of gravity originating from moving Jupiter. The light cone is a solution of the Maxwell equations 
and describes the null characteristics of a radio wave from a quasar. The gravitational deflection of light is seen at the point 
where the gravity cone of Jupiter passes through the observer. VLBI measures the Minkowski dot product $ = —kar" between 
two null vectors k°' and = a;" — Xj{s) at the point of observation directed to the observer from the quasar and from Jupiter, 
respectively. Had moving Jupiter not been detected at the retarded position on its world line, the speed of gravity would not 
be equal to the speed of light and the general theory of relativity would be violated. 
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FIG. 2: Light propagates from the quasar toward a VLBI station (observer) on the Earth. While the hght propagates, Jupiter 
is moving. General relativity predicts |3 that the light observed at the time t is deflected most strongly by Jupiter when it 
is located at the retarded position xj{sg) {sg = t — r/cg) for Cg = c, regardless of the direction of propagation of the light ray 
and the magnitude of the light-ray impact parameter with respect to Jupiter. This property of the gravitational deflection of 
light by a moving massive body can be used in order to measure the aberration of gravity force with respect to the aberration 
of light. The region of the sky which is imaged by the VLBI observations is smaller than the angle G, so that radio emission 
of Jupiter is not detected when the quasar is observed. 
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FIG. 3; The two null cones of the bi-metric theory are shown (see appendix). Retarded position of Jupiter due to the finite 
speed of gravity is on the gravity null-cone while Jupiter's radio position is on the light null-cone. The unperturbed direction 
to the quasar is defined by vector a°' lying on the light null-cone of the optical metric Qa/s- Gravity perturbs the light null- 
cone vector (j" and changes its direction to the gravity null-cone vector Gravitational field of moving Jupiter deflects the 
quasar's light from the retarded position defined with respect to observer by vector = a:" — Xj (s) which lies on the gravity 
null-cone. Perturbed eikonal ip ~ 2Mjln$, where $ = —lar°'. Components of the vector Z° are calculated from vector cr" 
and four-velocity of Jupiter u" which are known. The angle ^ is calculated assuming that the retarded position of Jupiter is 
defined by the JPL ephemerides and lies on the light null-cone. The observed value of the angle $ is fit to that defined by the 
JPL ephemerides by adjusting the parameter e defining the difference between the speed of light and gravity. This procedure 
measures the speed of gravity with respect to light within the accuracy of 20% for 10 pas astrometric tolerance in observing 
small angles in the sky with the phase-referencing VLBI technique. 



